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Introduction 

This report summarizes our investigations performed under NASA Grant NAG5-8058. 
The three-year research supported by the Geospace Sciences SR&T program (Iono- 
spheric, Thermospheric, and Mesospheric Physics) has been designed to investigate 
fluxes of energetic oxygen and nitrogen atoms in the terrestrial thermosphere. Fast 
atoms are produced due to absorption of the solar radiation and due to coupling be- 
tween the ionosphere and the neutral thermospheric gas. We have investigated the 
impact of hot oxygen and nitrogen atoms on the thermal balance, chemistry and radi- 
ation properties of the terrestrial thermosphere. Our calculations have been focused on 
the accurate quantitative description of the thermalization of O and N energetic atoms 
in collisions with atom and molecules of the ambient neutral gas. Upward fluxes of 
oxygen and nitrogen atoms, the rate of atmospheric heating by hot oxygen atoms, and 
the energy input into translational and rotational-vibrational degrees of atmospheric 
molecules have been evaluated. Altitude profiles of hot oxygen and nitrogen atoms have 
been analyzed and compared with available observational data. 

Energetic oxygen atoms in the terrestrial atmosphere have been investigated for 
decades, but insufficient information on the kinetics of fast atmospheric atoms has been 
a main obstacle for the interpretation of observational data and modeling of the hot 
geocorona. The recent development of accurate computational methods of the colli— 
sional kinetics is seen as an important step in the quantitative description of hot atoms 
in the thermosphere. Modeling of relaxation processes in the terrestrial atmosphere has 
incorporated data of recent observations, and theoretical predictions have been tested 
by new 7 laboratory measurements. 

The main results obtained under the present proposal are: 

1. The method of numerical solution of the kinetic equation, describing thermal- 
ization process, has been developed and has been used to compute non-Maxwellian 
distributions of atmospheric oxygen and nitrogen atoms. Fluxes of hot 0 and N atoms 
in the terrestrial thermosphere have been obtained in numerical solutions of the time- 
dependent Boltzmann equation. 
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2. The rates of energy-transfer collisions of energetic 0 and N atoms with ambient 
gas have been computed with accurate collisional cross sections. The angular and energy 
dependences of cross sections were obtained in quantum calculations. 

3. Thermalization times for energetic 0( 3 P), O^D), and N( 4 S) atoms and the atmo- 
spheric heating rate have been calculated at different altitudes for daytime and nighttime 
atmospheric conditions. 

4. Fractions of energetic oxygen 0( 3 P) atoms have been computed for altitudes of 
200-700 km at different solar conditions, and theoretical results have been compared 
with observations. 

5. Non-Maxwellian distributions of fast N( 4 S) atoms in the terrestrial atmosphere 
have been calculated with updated sources of nascent energetic nitrogen atoms. Com- 
puted distribution functions have been used for analysis of NO production rate in the 
terrestrial thermosphere. Theoretical results have been compared with observations. 

6. Parameters of the upward flux of energetic 0( 3 P) atoms have been determined in 
numerical calculations. The angular, energy, and altitude distributions of fast atoms in 
the upward flux near the exobase have been obtained from results of the Monte Carlo 
simulations of 0( 3 P)+0( 3 P) collisions carried out with realistic energy and angular 
dependent cross sections. 

7. The probabilities of escape from the thermosphere to the exosphere have been 
computed for ground state 0( 3 P) atoms produced at different altitudes below the 
exobase. 

The description of the theoretical methods and obtained results have been published 
in our articles and reported on AGU Meetings: 

Kharchenko V. and A. Dalgarno, Fluxes of fast atoms near the terrestrial thermo- 
sphere, Eos. Trans. AGU, 82(20), S293 (2001) 

Kharchenko V., A. Dalgarno, B. Zygelman, and J.-H. Yee, Energy transfer in col- 
lisions of oxygen atoms in the terrestrial atmosphere, J. Geophys. Res., 105, 24899 
( 2000 ) 

Kharchenko V., N. Balakrishnan, and A. Dalgarno, The upward flux of fast oxygen 
atoms in the terrestrial thermosphere, Eos. Trans. AGU, 81, S333 (2000) 

Balakrishnan N., E. Sergueeva, V. Kharchenko, and A. Dalgarno, Kinetics and ther- 
malization of hot N( 4 S) atoms in the upper atmosphere, J. Geophys. Res., 105, 18549 
( 2000 ) 

Kharchenko V., and A. Dalgarno, Fraction of energetic oxygen atoms in the upper 
terrestrial thermosphere, Eos. Trans. AGU, 81, F931 (2000) 

Balakrishnan N., V. Kharchenko, and A. Dalgarno, Translational energy relaxation 
of hot 0{ l D) atoms , J. Phys. Chem., 103, 3999 (1999) 

Balakrishnan N., V. Kharchenko, and A. Dalgarno, Translational energy relaxation 
of hot 0( l D) atoms in collisions with N 2 , Eos. Trans. AGU, 80, S239 (1999) 
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Results of our recent investigations of near exobase fluxes of energetic oxygen atoms 
and details of Monte Carlo simulations will be reported in 

Kharchenko V. and A. Dalgarno, Fast OfP) and O^D) atoms in the upper terres- 
trial thermosphere, J. Geophys. Res., to be submitted (2002). 


Differential and total cross sections for collisions of O 
and N atoms with the atmospheric gas 

Nascent energetic atoms in the terrestrial thermosphere are produced by various 
photo-chemical and ionic reactions, and atmospheric sources of hot O and N have been 
carefully studied [ Cotton et al. 1993', Gerard et al. 1995,' Richards et al. 1994, Hickey 
et al. 1995; Shematovich et al. 1992, 1994; P° x an< ^ A. Hac, 1997). Translational en- 
ergies of nascent fast atoms are larger on two - three orders than thermal energy of the 
ambient gas, and hot atoms slow down in multiple collisions with thermal atmospheric 
atoms and molecules - O, N 2 , and 02- Energy distributions of fast O and N atoms, the 
rate of the atmospheric heating by fast atoms, and the energy input into translational 
and rotational-vibrational degrees of the bath gas molecules depend strongly on colli- 
sional cross sections. Calculations of these cross sections are complicated due to quantal 
character of atom -atom and atom-molecule collisions at eV’s energies. Strong energy 
and angular dependences of realistic cross sections have to be taken into account in any 
evaluations of thermalization rates. 

Accurate calculations of cross sections 

Cross sections and thermalization rates for O + N 2 , N and N -I- N 2 , O 2 collisions have 
been obtained in our earlier calculations [ Kharchenko et al. 1997, 1998, Balakrishnan et 
al. 1998a, 1998b). Atom - molecule collisions are important in the low thermosphere, but 
collisions with atmospheric atomic oxygen dominate at altitudes above ~ 200 km. We 
performed accurate calculations of the energy transfer rate in 0( 3 P) +0( 3 P) collisions 
[Zygelman, Kharchenko, and Dalgarno 1999; Kharchenko et al. 2000], which are the 
main thermalization mechanism in the middle and upper thermosphere. Cross sections 
for the thermalization of fast metastable 0( ! D) atoms were also calculated [B alakrishnan 
et al. 1999] and compared with available experimental data [ Matsumi et al. 1994, 
Matsumi and Chowdhury 1996). 

W r e computed differential and total cross sections for elastic collisions of ground 
state oxygen atoms, using 18 electronic energy curves of 0 2 molecules [Zygelman et 
al. 1994a, b). Electronic states connected to the ground state of oxygen atoms are 
represented by singlet, triplet and quintet spin states. In Fig.l the computed cross 
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sections for the individual spin states of colliding oxygen atoms are shown as functions 
of the kinetic energy. 



ENERGY [eV] ENERGY [eV] 


Fig.l Fig-2 

Total and differential cross sections were calculated for all individual spin states, and 
results have been statistically averaged over all spins. The calculated statistical cross 
sections (solid curve) are well described by the Landau-Schiff approximation (dashed 
curve) as shown in Fig. 2 [ Landau and Lifshitz, 1991). The computed differential cross 
sections of 0 + 0 collisions are given in Fig. 3 as functions of the scattering angle for 
collision energies of 0.03 , 0.3 and 3 eV. 



Strong angular anisotropy of differential cr oss sections is induced by long-range van 
der Waals forces. The statistically averaged cross sections are similar to those for elastic 
collisions of excited oxygen atoms [Fee and Dalgarno 1985, 1987). 
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New values of the transport coefficients in oxygen gas 

Transport properties of the hot nitrogen and hydrogen atoms in the terrestrial atmo- 
sphere were calculated using accurate differential cross sections of elastic and inelastic 
collisions. Collisions of energetic oxygen and nitrogen atoms with thermal atomic oxy- 
gen are most important for the altitudes above 180 km. The reduced transport integrals, 
which are widely used for the simplified description of momentum and energy relaxation 
processes [Chapman and Cowling 1970], were calculated as functions of the gas temper- 
ature and were compared with previous theoretical and experimental data. In Fig-4 
results of our calculations of the a\ 1, 1] and a[ 2, 2] transport integrals in 0( P)+ 0( P) 
collisions are shown by solid curves, the dashed curves show the results of quasi-classical 
calculations by Levin et al. (1990), and the diamonds are the results of classical trajec- 
tory calculations by Holland and Biolsi (1988). 

The self-diffusion coefficient D s and the viscosity 77 [Chapman and Cowling 1970; 
Ferziger and Kaper 1972] have been calculated for the temperature interval of 100 K - 
30000 K and have been compared with available data [Dalgarno and Smith 1962, Levin 
et al. 1990]. Calculated transport coefficients D s and 77 may be fitted by interpolation 
formulas 

n ~ - — * T 0 - 8 C7n 2 s' 1 , 77 — 2.72 * T° 76 * 10~ 6 g cm' 1 *- 1 , 

n 0 

where T and n 0 are the temperature and the density of the atmospheric oxygen gas. 


Energy relaxation of hot O and N atoms in the ter- 
restrial thermosphere 

Accurate differential cross sections of O and N collisions with atmospheric atoms 
and molecules have been used in evaluations of the energy relaxation kernel B(E f \Ei), 
which gives the rate of transitions from the state with initial energy E x to the final state 
with translational energy E f [Oppenheim et al. 197l\. The collision kernel has been 
calculated using analytical and numerical methods, developed in our investigations of 
the thermalization of hot nitrogen and oxygen atoms [ Kharchenko et al. 199 7, 1998, 
2000 ; Balaknshnan et al. 1998a, 1999, 2000]. In Fig.5 the calculated rates of E { E f 
transitions at different initial energies of 0.1, 1, 2.5 and 4.5 eY are shown for oxygen 
atoms thermalizing in the oxygen gas with the temperature of 850 K and the density of 
3*10 9 cm' 3 [Kharchenko et al. 2000]. The energy transition rate for the atmospheric gas 
mixture is calculated as a statistically weighted sum of the collisional kernels B l {Ef\E i ) 
of the individual atmospheric gases (i= O, N 2 , and 0 2 )- The rate of energj transitions 
B(Ef\Ei) strongly depends on the altitude, because of the altitude dependence of the gas 
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density and composition. We have calculated the energy-transfer rates for fast oxygen 
atoms in the terrestrial atmosphere between 100 and 600 km. The rate of atmospheric 
heating normalized to a single hot O atom with given initial energy E has been also 
computed at different altitudes. Results are shown in Fig. 6 as function of kinetic energy 
E for the altitude range from 100 km to 400 km. 




Fig-5 


Fig. 6 


Times of thermalization of fast atoms depend on their initial energies. Distribution 
of the kinetic energies of nascent fast atoms produced by the various photo-chemical 
and ionic reactions [ Richards et al. 1994 ', Hickey et al. 1995; Fox and Hac 199 7] are 
shown in Fig.7 by dashed curves at 200, 300, and 400 km. In the same figure, the 
thermalization times of energetic O atoms at different altitudes were calculated with 
accurate quantal cross sections [Kharchenko et al. 2000}. 

We calculated the numbers of collisions, w'hich are required for the thermalization 
of energetic O and N atoms in the atmospheric gas at different altitudes. The number 
of thermalizing collisions increases dramatically, if the initial energy of the nascent hot 
atom increases. Because of the strong anisotropy of differential cross sections the hot 
atoms lose their energy mostly by small amounts in long distance collisions with ambient 
atoms and molecules. 

Energy relaxation process involves dozens of collisions, and the thermalization time 
increases significantly with the altitude. In Fig. / the thermalization times for the atoms 
with different initial energies •E’feV'] are shown by solid lines for altitudes 200-400 km. 
Thermalization times of O and N atoms and rates of their energy losses were calculated 
for daytime and nighttime conditions of the terrestrial thermosphere. Rates of energy 
losses of fast 0( 3 P) and N( 4 S) in collisions with atmospheric atoms and molecules were 
calculated for gas temperatures, densities, and compositions corresponding to the pa- 
rameters of the terrestrial thermosphere between 100 - 600 km. 
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Fig. 7 


Fig-8 


Secondary hot oxygen atoms in the upper thermosphere 

In the middle and high regions of the terrestrial thermosphere, where atomic oxygen 
is the main constituent of the atmospheric gas, secondary hot 0( 3 P) atoms arise in recoil 
collisions of nascent energetic atoms. Due to high efficiency of the energy transfer in O - 
O collisions the recoil atoms may receive significant fraction of the kinetic energy, which 
was carried initially by nascent fast atoms. Recoil oxygen atoms play an important 
role in the formation of upward fluxes of energetic 0 atoms. We calculated the energy 
distribution of the secondary energetic oxygen atoms in binary collisions with fast N 
and O. These data have been used in our evaluations of the kinetics of hot atoms and 
in computing of atomic upward fluxes near the exobase. 

Energy spectra of the resulting secondary energetic atoms have been obtained in 
quantum mechanical calculations of the energy transfer in O + O collisions [ Kharchenko 
et al. 200]. In Fig.8 the spectra of recoil oxygen atoms computed with accurate energy- 
transfer cross sections of O + O collisions are shown for different initial energies E t of 
nascent hot O atoms. Distribution functions of recoil atoms are normalized to unity. 
Number of recoil atoms increases exponentially with time, but their spectra converge 
asymptotically to the thermal distribution. 

N on-Maxwellian distributions of fast O and N atoms 

Calculations of non-Maxwellian distribution functions of fast O and N atoms are 
based on the solution of the Boltzmann equation describing the thermalization of the 
nascent energetic atoms [ Shizgal and Blackmore 1986; Shizgal 198 7; Stamnes et al. 
1991; Lie-Svendsen et al. 1991, 1992; Shematovich et al. 1992, 1994, 1999; Bisikalo et 
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al. 1995; Sharma et al. 1996 ; Gerard et al. 1995, 1997; Kharchenko et al. 1997, 1998, 
2000; Balakrishnan et al. 1998a,b, 1999, 2000; Hubert et al. 1999; Kharchenko and 
Dalgarno 2000a, b). Solutions of the Boltzmann equation or Monte Carlo simulations of 
the energy relaxation process provide accurate non-Maxwellian distribution functions 
and escape fluxes of energetic atoms, if realistic differential cross sections of thermaliz- 
ing collisions are used. Reported in the literature solutions of the kinetic equation and 
Monte Carlo simulations were mostly carried out for the hard sphere model, which as- 
sumes collisional cross sections independent of the energy and scattering angles [Shizgal 
and Blackmore 1986; Shizgal 1987; Stamnes et al. 1991; Lie-Svendsen et al. 1991, 1992; 
Shematovich et al. 1992, 1994, 1999; Sharma et al. 1996; Marov et al. 1997; Hubert et 
al. 1999}. With this model the calculation procedure is simplified significantly. How- 
ever, the hard sphere approximation does not describe adequately the energy -transfer 
collisions, because realistic cross sections of atom-atom and atom-molecule collisions are 
stronglv focused in the forward direction and depend on the collision energy [ Kozyra et 
al. 1982; Yee and Dalgarno 1985, 1987; Gerard et al. 1997; Kharchenko et al. 1997, 
1998, 2000; Balakrishnan et al. 1998a,b, 1999, 2000]. 

The new approach, which incorporates the quantum treatment of elastic and inelastic 
collisions with the solution of the Boltzmann equations, was developed for calculations 
of the energy distribution function of fast nitrogen atoms in the terrestrial thermosphere 
[Kharchenko et al. 1997, 1998 ; Balakrishnan et al. 1998a, 2000]. Recently it has been 
used to compute the energy relaxation of the ground-state and excited 0 and N atoms 
in the atmosphere [ Balakrishnan et gl, 1998b, 1999; Zygelman et al. 1999; Kharchenko 
et al. 2000; Balakrishnan et al. 2000; Kharchenko and Dalgarno 2000b}. 


The kinetic equations for thermalization of fast O and N atoms 

Thermalization of hot atoms in the terrestrial thermosphere is governed by elastic 
and inelastic collisions with the thermal 0( 3 P) atoms and N 2 and 0 2 molecules. The 
computing scheme for calculations of the energy relaxation and the upward fluxes of 
fast O and N atoms are complicated by the presence of nascent metastable 0 and N 
atoms. Quenching collisions of excited atoms with atmospheric gas generate fast ground 
state O and N atoms, and distributions of ground state and excited atoms should be 
calculated self-consistently. Solutions of the Boltzmann equations allow to calculate up- 
ward fluxes of metastable and ground states atoms in consistent conditions. From the 
steady' state solutions the altitude profiles and energy distributions of thermospheric 
oxygen and nitrogen were determined and compared with available observations and 
theoretical predictions [Cotton et al. 1993 ; Gerard et al. 1995, 1997; Oliver 1997, 
Swaminathan et al. 1998; Hubert et, al. 1999; Oliver and Schoendorf 1999; Litvin et al. 
2000; Kharchenko et al. 2000; Kharchenko and Dalgarno 2000b, 2001; Balakrishnan et 
al. 2000; Schoendorf et al. 2000]. 
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Collisions between fast particles may be neglected in investigations of upward fluxes 
of energetic atoms, because hot atoms are only a minor fraction below 700-800 km 
[i Cotton et al. 1993; Bisikalo et al. 1995; Gerard et al. 1995; Shematovich et al. 
1999; Kharchenko and Dalgarno 2000b, 2001}. At higher altitudes atomic fluxes are 
well described by collisionless transport. The system of the Boltzmann equations for 
the momentum distribution functions F,(p,z,t) of hot atoms with the identities of 
i = {0{ 3 P),0{ l D), Oi'S), N( 4 S), N( 2 D)} are given by 

z, t) + m t g * VpFi(p , z, t) + 2 > *) = 

J dpxFi(pi,z,t)Wi(p,pi,n z ) - oji(p,z)F t (p,z,t) + Qi{p,z,t), 

where m* is the mass of the fast atom, g is the free-fall acceleration, z is the vertical axes, 
Q t (p,z,t) is the intensity of the source reactions, w t (p,z) is the frequency of elastic and 
inelastic collisions between fast atoms and the ambient thermal gas, and Wi(p,pi,n z ) is 
the momentum relaxation rate depending on the differential cross sections of elastic and 
inelastic collisions and on the composition, temperature and density n z = n(z) of the 
atmospheric gas. The collisions, describing quenching of metastable O and N atoms, are 
taken into account in the source functions of fast 0( 3 P) and N( 4 S) atoms [Gerard et al. 
1995; Hickey et al. 1995; Shematovich et al. 1999; Kharchenko and Dalgarno 2000b}. 
In the upper thermosphere, secondary energetic recoil 0( 3 P) atoms are included in the 
source function with the spectra shown in Fig. 8 [ Zygelman et al. 1999; Kharchenko 
et al. 2000}. Analytical expressions of the collision kernel via differential and doubly 
differential cross sections of elastic and inelastic collisions, which have been used in our 
numerical computations, were reported in articles [Kharchenko et al. 1997, 1998, Bal- 
akrishnan et al. 1998a,b; Shizgal 1999}. 

An alternative method for the numerical solution of the kinetic equation involves 
Monte Carlo simulations of multiple collisions between energetic atoms and the thermal 
atmospheric gas [Shematovich et al. 1994, 1999; Bisikalo et al. 1995; Gerard et al. 
1995; Hubert et al. 1999; Kharchenko and Dalgarno 2001}. 


Monte Carlo simulations of upward traveling fluxes 

Monte Carlo simulations may be a very efficient method of calculations of upward 
fluxes of fast atoms, once accurate energy and angular dependences of scattering cross 
sections have been used for the description of stochastic binary collisions. Monte Carlo 
simulations were carried out in our recent calculations to determine the energy and an- 
gular distributions in the vertical flux of hot 0( 3 P) atoms [Kharchenko and Dalgarno 
2001 }. 

• An ensemble of hot oxygen atoms with energies E{ n has been generated at different 
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initial altitudes H in . Initial velocity directions of nascent fast atoms have been given by 
random numbers, and initial energy distributions have been compiled form the known 

spectra of source reactions. - - ----- 

• Collisions of ground state hot atoms with the thermal gas were described by four 
stochastic variables: j - the index of the identity of the bath gas atoms or molecules 
(O or N 2 for the atmosphere above 200 km ); R - the projectile radius- vector with the 
altitude projection of z; 6 - the angle between the velocity V and the vertical axes z\ x 
- the Center of the Mass Frame (CMF) scattering angle (0 < x < *0; $ - the azimuthal 
angle between the scattering and [V, e z \ planes (0 < $ < 2tt). 

• The altitude, angular, and energy distributions of hot oxygen atoms have been ob- 
tained from the results of Monte Carlo simulations of multiple collisions between fast 
atoms and atoms of the ambient gas. 

These distributions have been used to compute the fractions of energetic atoms escaping 
from the thermosphere to the exosphere [Lie-Svendsen et al. 1992; Hubert et al. 1999, 
Kharchenko and Dalgamo 2001}. In Fig.9 the calculated probabilities to escape into the 
exosphere (above 700 km) are shown as function of the initial altitude, where nascent 
fast 0( 3 P) atoms have been produced. For atoms with the initial energy of 5 eV initial 
altitudes have been taken from 400 km to 700 km. 


escape from the thermosphere to the exosphere 

0( 3 P) hot atoms , H M1 = 700km , N aloms = 10 



Calculations have been carried out for the quantally computed anisotropic cross 
sections and for the hard sphere model. Monte Carlo simulations have been used for 
computing upward fluxes of energetic O and N atoms in the upper thermosphere. 
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Time-dependent calculations of the energy relaxation. 

The time-dependent energy distributions of nascent hot atoms provide detailed in- 
formation on the relaxation rate and local heating of the ambient gas in conditions of 
variable sources of fast atoms. Time-dependent thermalization of energetic O and 
atoms in the exobase region have been studied at variable magnetospheric and iono- 
spheric conditions, and at different levels of solar activity. For the solution of the 
time-dependent Boltzmann equation we developed an efficient numerical method based 
on the matrix representation of the evolution operator via the energy relaxation ker- 
nel [Zygelman et al. 1999; Kharchenko and Dalgarno 2000b, 2001}. The results of our 
investigations of the energy relaxation of metastable and ground state O atoms and 
ground state N atoms have been reported recently [ Balakrishnan et al. 1999; Zygelman 
et al. 1999; Kharchenko and Dalgarno 2000a, b], and the theoretical data have been 
compared with the results of laboratory measurements [ Matsumi et al. 1994; Matsumi 
and Chowdhury 1996; Balakrishnan et al. 1999; Caledonia 2000]. Numerical solutions 
of the kinetic equation with realistic differential cross sections provide reliable informa- 
tion on the energy relaxation. In Fig.10 the results of our calculations of the energy loss 
of O^D) atoms in the N 2 atmosphere are shown to be in excellent agreement with the 
experimental data [ Matsumi et al. 1994; Matsumi and Chowdhury 1996; Balakrishnan 
ct al . 1999]. 

Time-dependent solutions of the Boltzmann equation have been used for investiga- 
tions of the time evolution of energy distributions of 0( 3 P) and N( 4 S) atoms in the 
terrestrial thermosphere [Zygelman et al. 1999; Kharchenko et al. 2000; Kharchenko 
and Dalgarno, 2000a, b}. Fig. 11 shows the time-dependent energy ^distributions of ener- 
getic oxygen atoms in the terrestrial thermosphere at 0.03 s, 0.17 s, 0.34 s and 1.35 s 
after they have been produced. All distributions are normalized to unity 




Fig.10 Fig. 11 

We included in our calculations more then twenty sources of nascent hot oxygen 
atoms [ Shematovich et al. 1994; Gerard et al. 1995; Hickey et al. 1995 ]. 
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Computational method, which has been developed in the present proposal, could 
be used in calculations of energy distributions of metastable and excited atoms. The 
energy and velocity relaxation of nascent metastable atoms may be determined fiom 
solutions of the time-dependent Boltzmann equations, and resulting distributions could 
be used in calculations of the Doppler profiles of OpD/S) emissions detected in the 
atmospheric observations [Schmitt et al. 1981, 1982; Killeen and Hayes 1981, 1983, 
1984; Shematovich et al. 1999]. 
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